Coulomb drag between helical edge states 
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FIG. 1. (color online) Schematic of a drag measurement be- 
tween two QSH systems. A current ii is driven along the 
upper edge of the lower QSH system and through electron- 
electron interactions a voltage V2 is induced in the lower edge 
of the upper QSH system. A magnetic field h is applied in 
the plane of wires, perpendicular to the spin quantization axis 
(assumed perpendicular to the plane of QSH systems). Time- 
reversed Kramer's pairs are indicated for the two edges. A 
QSH on top of QSH geometry could also be used. 
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We theoretically investigate the Coulomb drag between the edge states of two quantum spin Hall 
systems. Using an interacting theory of the one-dimensional helical edge modes, we show that 
the drag vanishes at second order in the inter-edge interaction, where it is typically finite in other 
systems, due to the absence of backscattering within the edges. However, in the presence of a small 
external magnetic field the drag is finite and scales as the fourth power of the magnetic field, a 
behavior that sharply distinguishes it from other systems. We obtain the temperature dependence 
of the drag for regimes of both linear and quadratic edge dispersion in the presence of a finite field. 



Introduction-Topologies! phases of matter have at- 
tracted interest because their quantum properties are 
robust to many material imperfections ji — In particu- 
lar, the quantum spin Hall (QSH) system has recently 
emerged 4 ^— as a time-reversal invariant counterpart to 
the integer quantum Hall effector— Shortly after the 
prediction 1 ^ that the HgTe/ (Hg,Cd)Te quantum well sys- 
tem should exhibit a QSH state, the experimental obser- 
vation was made 1 - 1 - and further confirmation followed. 12 

The QSH state has an insulating bulk and metallic 
edge states composed of an odd number of Kramer's pairs 
of electrons. A Z2 invariant distinguishes the topolog- 
ical insulators with time-reversal symmetry from their 
"trivial" counterparts^— The simplest topologically non- 
trivial case is a single Kramer's pair on the edge. Due 
to the spin-orbit coupling that drives the QSH state, 
the spin of an electron on the edge is correlated with 
its momentum. This property leads to an absence of 
back-scattering from weak non-magnetic impurities and 
therefore prevents Anderson localization on the edge of 
the QSH system,! 3 *! 4 . 

The gapless edge modes of the QSH system are com- 
monly referred to as a helical liquid (HL)J^ The stability 
of the HL to interactions ; 13 ! 14 and magnetic disorde r 13 ! 15 
has been investigated, as has its response to "pinch- 
ing" the sample into a point contact**— or related 
geometries^ Properties of superconducting-QSH hybrid 
structures were investigated as well<^ When two HL (of 
different QSH systems) are allowed to interact with each 
other, a novel one-dimensional correlated state is formed 
at the lowest energies^ 

In this paper we study the Coulomb drag between two 
QSH systems, each having one Kramer's pair as shown 
schematically in FigJT] The drag experiment we dis- 
cuss should be carried out at energy (temperature) scales 
above which the inter-edge correlated state forms; 21 - but 
below the bulk energy gap of the QSH system. In a 
Coulomb drag experiment current is driven in an "ac- 
tive" wire/edge and voltage is measured in a "passive" 
wire/edge. The Coulomb interaction between electrons 
in different system results in a momentum transfer be- 
tween the two systems and produces a voltage drop in 
the "passive" system. The drag is characterized by the 



drag resistivity, 
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where V% is voltage induced in the "passive" system by 
the current 1\ driven in "active" system. Here e is the 
electron charge, h is a Planck's constant, and L is the 
length of the edge along which momentum is transferred. 

Coulomb drag in non-HL one-dimensional sys- 
tems has been studied both theoretically 2 ^— and 
experimentally^ - — The HL can be viewed as a spin- 
less Luttinger liquid (because it has the same number 
of degrees of freedom) without backscattering; since it 
is known that the backscattering governs the drag be- 
tween systems with linear dispersions ; 22 ' 25 ! 27 one can not 
expect drag between helical liquids. However, in this pa- 
per we show that an applied Zecman field h opens up a 
backscattering process in the HL, and results m rr> oc h 4 
for a linear spectrum. We also compute the temperature 
dependence of the drag over a range of temperature and 
field values. 

The model-In second order perturbation theory in the 
interwire Coulomb interaction, the Coulomb drag is given 
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by^ 
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where SsTIi(a;, q) is the imaginary part of the retarded 
density-density correlation function of wire i — 1(2), rii 
is electron density of wire i, T is the temperature, and 
U\2{q) is the Fourier transform of the interwire Coulomb 
interaction which is cut off at short distances by the in- 
terwire separation d. 

We consider two identical QSH systems, each with one 
Kramer's pair on its edge, as shown in FigfTJ As we noted 
earlier, if the spectrum is linear there is no contribution 
to the drag from forward scattering, and back scattering 
is forbidden by time-reversal symmetry. Therefore, one 
must break time-reversal symmetry in order to open up 
a backward scattering channel (unless there are magnetic 
impurities present) for a generic Dirac edge mode. Our 
Hamiltonian for a single HL in the presence of a Zeeman 
field is 

H = J dx${x) (vp x a z + ha x - p) ip(x), (3) 

where v is the edge velocity, p x = —id x , p, is the Fermi 
energy (which can be adjusted by gating the system), 
and a ZiX are Pauli spin matricies describing the spin de- 
gree of freedom. A Zeeman field h pointing in the re- 
direction opens up a gap in Dirac spectrum and tilts 
the spins away from the z-axis. The edge dispersion is 
e± = ± sj v 2 p 2 + h 2 — p. We assume that Fermi energy 
is in the upper band (// > 0) so that the properties of 
system are determined by the e + band over the energy 
scales of interest. The wavefunction of electrons in the 
e + band is 



ev* / cos( 7p /2) + sm( 7p /2 

V>+W = cos( 7p ) 

V2 \ (cos( 7p /2)+ S in( 7p /2)) 

where t„ = arctan(^). We studv (T- 



e ipx U p , (4) 



where 7 P = arctan(^). We study ((4]) in the limit of 
large p, (small h) where the spectrum can be approxi- 
mated as linear, and also in the opposite limit where the 
spectrum is approximately quadratic (i.e., p close to the 
band "bottom"). See Figs, [land El 

Regime of linear spectrum-We first consider the case 
(i — h > h, and linearize the spectrum near the Fermi 
energy, e + = v\p x \ — /i (see Fig. [5]), in order to use 
standard bosonization procedures.— We express the elec- 
tron operator as a sum of left- and right- moving states: 
%l> + = tPr{x) + i/jl( x ), where R{L) stands for right (left) 
movers. The non- interacting Hamiltonian can then be 
written: 



dx 



Wr{x)p+^r{x) + ipl(x)p-ipL(x) 



(5) 



where p± = ±vp x — /i. We assume intrawire interactions 
have the form 

Hmt = 2 dx dx' U(x - x')p(x)p(x'), (6) 




FIG. 2. (color online) Temperature dependence of the drag 
in the regime of small h where the spectrum may be approx- 
imated as linear, as shown in the inset. T* is the tempera- 
ture at which the wires begin to "lock" to each other.— For 
T > T*, we find r D oc h 4 T 4K ~ 3 where K is the Luttinger 
parameter in the charge sector. 



where U(x-x') is the intrawire Coulomb interaction, and 
p{x) is the electron density: 



p(x) = ^ R i! R + tplip L + cos(7 P ) (V^l + V'lV'ij ) , (7) 



which contains cross terms due to the presence of the 
magnetic field. In terms of bosonic fields <fi( x ) ancl @( x ), 
if)R and are expressed as^ 2 . 
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where ??_r(l) are Klein factors, and a is a short-distance 
cut-off. The electron density in terms of bosonic fields 
takes the form 



p(x) 
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d x (j)(x) 



■sm(2p F x- 2cj){x)). (10) 



Substituting this expression into (J6j) , we find 

las i i(.)-^,W») [ dxmix)r , (11 ) 
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where U (0) and U (2pf) are the zero and 2pF momentum 
parts of the interaction, respectively. Note that the 2pF 
part has a cos 2 (7p F ) factor, which is proportional to h 2 
for small h. The full Hamiltonian then becomes 



H 



1 

2^ 



dx \v(d x e(x)) 2 + (v + g) (d x 0(x)) 2 ] , (12) 



where g = (£7(0) — cos 2 (7 Pf )E/(2pf))/7t. We observe 
that the Hamiltonian of an interacting HL in a Zeeman 
field is equivalent to a spinless Luttinger liquid where the 
strength of backscattering depends on the Zeeman field. 
Similar results were obtained in studies of Luttinger liq- 
uids with Rashba spin-orbit coupling and a Zeeman mag- 
netic field^^ 

We now give an expression for the retarded density- 
density correlation function. Because the backscattering 
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governs the Coulomb drag when the dispersion is linear 
(as it is in a Luttinger liquid), we only need the 2p F 
part of the retarded density-density correlation function. 
Since our model has reduced to a spinless Luttinger liq- 
uid, the calculation is standard^ 3 - 2 - 
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Ii(q + 2p F ,Uj) + Ii(q-2pF,Uj) , (13) 



with f[(g, oj) is given by 
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where F(q,uj) = B(-i^(uj-uq) + f,l-K)B(-i^{ 
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FIG. 3. (color online) Temperature dependence of the drag 
in the regime of small where the spectrum is approximately 

2 

quadratic, as shown in the inset, and < n — h < . Note 
the non-monotonic temperature dependence^ above T* . For 
the dependence of ru on the Zeeman field h in each region of 

to 



l — K), [5 = 1/T, u = \Jv(v + g), K = v/u, and temperature, see the text. The second crossover from T 



B{x, y) is the Beta function. The parameter D is 



D = cos 2 (j PF ) sin(7rif ) 
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T' 1 occurs for T 
wires. 



where d is the distance between 



With (|T4|) in hand, the drag resistivity is readily com- 
puted from @: ^SF(q± 2p F ,u)) is sharply peaked about 
q = I f2p F with peak widths proportional to tempera- 
ture. Since the momentum integration in ([2} runs from 
to oo, we neglect the H(q + 2p F ,uj) contribution. Taking 
the imaginary part of the retarded density-density corre- 
lation function and assuming identical helical liquids we 
obtain 
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of states, n — X/irv. Extracting the temperature and 
magnetic field dependence using (fT5")) , we find 



r D oc 



h A T 4K - 3 . 
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Eq. (IT7)l is one of the central results of the paper. This re- 
sult is valid at temperatures larger than T* , below which 
the drag begins to exhibit an exponential dependence on 

temperature^^ Since T* ~ ^e" 1 ^ depends 2 ^ on the 
backscattering via K, it will also depend on the Zeeman 
field via the dependence of K on h. 

By contrast, in a spinful Luttinger liquid the mag- 
netic field only enters the interaction constant in the spin 
channel and therefore the drag is only (weakly) depen- 
dent on magnetic field through the interaction parameter 
appearing in an exponent to the temperature. There- 
fore, Coulomb drag can be used as a method for experi- 
mental verification of the HL, complementing the earlier 
studies . 13 ' 15 " — We note that a spin-Coulomb drag effect 
in which two density mismatched Luttinger liquids can 
be brought into more favorable kinematic conditions for 
enhanced drag effects has also been studied. 35 To com- 
plete our analysis of drag between two HL, we turn to 
the case when the spectrum is approximately quadratic. 



Regime oj quadratic spectrum-Whcn < fi — h < h, 
the spectrum of upper band is approximately e + = 

(<" — 'O' as seen m El In this section, 

we simplify the problem by assuming no intrawire elec- 
tron interactions. The neglect of weak interactions in the 
regime of a quadratic dispersion has been shown to have 
no effect on the temperature dependence of the dragi^ 5 . 
The imaginary part of the retarded density-density cor- 
relation function is 



4v 2 q 
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cosh(^^)cosh(- 



-(Po+g) 



2T 



) 

(18) 



up* 



which 



where p Q = -\q + 4^, and f+(p,q) = u^u p+q , 
we assume to be approximately equal to one. One also 
needs to take into account restrictions on u> defined by 
e + (po) < and e + (p + q) > which will give 
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h) > — > + yj2h{u-h). (19) 
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Plugging (fTgj) into ^ and evaluating the integrals we 



obtain the following results. When /i — h > 



2<r 2 u 



[here 



d is the inter-edge separation distance and 1/d serves as 
high momentum cut-off to U\2{q) in the q integration in 

©] and for small temperatures T < fg\ r^T' 



h 5 



rD 



(20) 



while at large temperatures T > \J~^jf, 
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When ji — h < and for small temperatures T < 
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while at large temperatures T > % 
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Here g 1 = —7 + ln(2) (7 w 0.5772 is Euler's constant) 
is an estimate of interwire Coulomb interaction at small 
momentum. The density of states, n = — , 1 , . The re- 

suits are summarized in Fig. [3] We emphasize that in ob- 
taining these results we have not considered effects of in- 
terband (intra-edge) particle-hole excitations. These ex- 
citations will result in Fermi edge singularity physics . 36 ! 37 
Summary and Discussion-We studied the Coulomb 
drag between identical one dimensional helical liquids. 
We showed that the helical liquid can be mapped to 
a spinless Luttinger liquid where backscattering is pro- 
hibited. Since backscattering governs the drag between 
one-dimensional liquids with linear dispersion, there is no 
Coulomb drag unless there is a nonlinearity in the spec- 
trum. Nonlinearity in the spectrum gives rise to a small 
momentum scattering contribution to the drag which has 
a T 2 temperature dependence^ Our calculations con- 
firm these results for a nonlinear spectrum. 



For a linear spectrum, the application of a Zeeman field 
opens up backscattering processes which are proportional 
to the square of the magnetic field at small fields, see (TPO)) 
and ([TT]) . We also showed that when the magnetic field 
is small (/1 large) and the spectrum can be approximated 
as linear, the Hamiltonian of a helical liquid with a mag- 
netic field is identical to that of a spinless Luttinger liq- 
uid with a magnetic field dependent backscattering term 
(IT2l . In this case, the Coulomb drag becomes propor- 
tional to the fourth power of magnetic field p7|) which 
is distinct from the case of Luttinger liquids where the 
magnetic field enters only via the interaction constant in 
the spin channel. 

For completeness we studied the case when the spec- 
trum of a helical liquid in a magnetic field can not 
be approximated as linear, but is rather approximately 
quadratic (valid for a strong magnetic field) . Expressions 
for the Coulomb drag in this case are given by (|2T))) - (|2"3"1) . 
Finally, we note that the presence of a few magnetic im- 
purities on the edge of a QSH system would allow a finite 
drag contribution even in the absence of applied magnetic 
fields since they would allow backscattering. Inclusion of 
Rashba coupling would not affect our results, provided 
the zero-field case is still adiabatically connected to the 
topologically non-trivial state. 
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